MONOTONE JACOBI PARAMETERS AND 
NON-SZEGO WEIGHTS 



YURY KREIMERi, YORAM LAST^ ^ AND BARRY SIMONY'S 



Abstract. We relate asymptotics of Jacobi parameters to asymp- 
totics of the spectral weights near the edges. Typical of our re- 
sults is that for a„ = 1, 5„ = —Cn~^ (0 < /? < |), one has 
dfj,{x) = w{x)dx on (—2,2), and near x = 2, w{x) — e^^'^'^^^ 
where 




1. Introduction 



Since the earliest days of the general theory of orthogonal polyno- 
mials on the real line (OPRL), it has been known that a key role is 
played by the Szego condition [38] that if 



where w is supported on [—2, 2] (we follow the spectral theorists' con- 
vention related to a„ — > 1, 6„ ^ rather than the [—1, 1] tradition in 
the OP literature), then 



In this paper, we will examine asymptotics of log{w{x)) for typical 
cases where (11. 2p fails. Recall [39l El El [3ll [3l] that, given /x, one 
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dfi{x) = w{x) dx + dfis 



(1.1) 




(1.2) 
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can define monic orthogonal and orthonormal polynomials Pn{x,dfi), 
Pn{x,dfi) and Jacobi parameters {a„,6„}^]^ by (6„ real, a„ > 0) 

XPn{x) = an+lPn+l{x) + hn+lPn{x) + anPn-l{x) (1.3) 

and 

||P„|| =ai---an (1.4) 
Favard's theorem (see, e.g., [311 [S^) asserts a one-one correspondence 
between /i's of compact but infinite support and bounded sets of a^'s 
and 6„'s. Moreover, by Weyl's theorem, if a„ — 1, 6„ — > 0, then the 
essential support of d/i is [—2,2]. 

Roughly speaking, the boundary for (11.21) to hold is a„ — 1, 6„ de- 
caying faster than 0{n~^). Explicitly, Killip and Simon [TT] proved 
a conjecture of Nevai [21] that Yl^=ii\'^n ~ 1| + \bn\) < oo ^ (11. 2p . 
and there are examples of Pollaczek |25l [261 [2Z| where (II. 2p fails be- 
cause \og{w{x)) ~ (4 — x^)~2 near a; = ±2 and 6„ = 0, a„ = 
1 -Cn-i + 0(n-2). 

Killip-Simon [TT] discovered a relevant weaker condition than (11.21) 
they called the quasi-Szego condition: 

y log(w(x))(4 — x^)5 c/x > — oo (1.5) 

and they proved that 

oo 

ffT3D+ (|a;| -2)t < oo ^ ^|a„-lp + |6„P < oo (1.6) 

3;Gsupp(/^)\[-2,2] n=l 

Our cases will include situations where (II. 5p and (11.60 fail. 

It is known (see [IHllSOl [2ll [22l [29l|40]) that when Y.n=Mn - Ip + 
= oo, djji can stop having an a.c. component, so we will need an 
additional condition. What we will use is 

Theorem 1.1. //a„ 1, hn — > 0, and 

oo 

^ |a„+i - a„| + - 6^1 < oo (1.7) 

n=l 

then (II. ip holds where w{x) is continuous on (—2, 2) and strictly posi- 
tive there. Moreover, dfig is supported on M \ (—2,2). 

The continuum Schrodinger analog of this is a theorem of Weidmann 
[H]; for OPRL, it is due to Dombrowski-Nevai [4j (see also [I2l[8l[32]). 
Most references do not discuss continuity of w but it holds; for example, 
it follows immediately from Theorem 1 of [3], since w can be obtained 
as a uniform limit of continuous functions on any closed subinterval of 
(-2,2). 
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In fact, we will focus on cases where {cin} and are monotone, 
so (11.71) is automatic. Typical is 

an = l hn = -Cn-" (1.8) 

where, roughly speaking, we will prove w{x) is singular at a; = 2 (i.e., 
the integral in (11. 5p diverges there) with 

w{x) = e-2«(-) (1.9) 
Q(s) ~Ci(2-x)5-^ (1.10) 

Indeed, in Section [Sj we will obtain for ( 11.81) an asymptotic series for 
Q[x) near x = 2 up to terms of 0(log(2 — x)); see ( 15.321) . 

Our interest in these problems was stimulated by a recent paper 
of Levin-Lubinsky [18] and their related earlier works on non-Szego 
weights [ini [T7]. They study the problem inverse to ours, namely, 
going from w (or Q) to a^, hn (which they call Bn). Unfortunately, 
they do not obtain even leading order asymptotics for an.hn if Q{x) 
has the form (ll.lOp but instead require 

Q(a;)~exp,(l-a:2)-" (1.11) 

with exp^(x) = exp(expj;,_i(x)) and exp]^(x) = e^. We will obtain 
inverse results to theirs in Section [51 We note that [16j does have 
asymptotics on the Rakhmanov-Mhaskar-Saff numbers when (11.101) 
holds and that their asymptotics should be connected to asymptotics 
of a„, bn- 

It is hard to imagine strict if and only if results on Q{x) to a„, hn 
since there will typically be side conditions {an,bn monotone and/or 
convex in n or Q{x) convex) that may not strictly carry over, but it 
is comforting (even with side conditions) to get results in both direc- 
tions. It would be interesting to show that (11.91) and (ll.lOp (with extra 
conditions) lead to estimates on a„,6„ with |a„ — 1| + \bn\ = 0{n~^). 
We suspect, with analyticity assumptions on Q, that this might be 
accessible with Riemann-Hilbert techniques. 

Our key to going from (a„, 6„) to {w, Q) is Carmona's formula that 
relates dfi to the growth of p„(a;), namely. 

Theorem 1.2. If Pn are the orthonormal polynomials for a measure 
djji, then du^'^^ — ^ dfx where 

^"'"^^'^^^-(a2p„(x)^ + pL,(x)) ^^-^2) 

The continuum analog of this result is due to Carmona [1]. This 
theorem when a„ = 1 is stated without proof in Last-Simon [T^ and 
later (with proof) in Krutikov-Remling [13\ and Simon [33] . It implies: 
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Corollary 1.3. Suppose uniformly on some interval [a, (3], we have for 
strictly positive continuous functions f± (x) that 

7r"V-(a;) < \imM{alpn{xy + Pn~i{x)^) 

< limsup(a„p„(x) ) < vr" 

Then dfi is purely absolutely continuous on {a,f3) and 

^ < w{x) < (1.14) 



Uix) - ^ ^ - /„(x) 
there. In particular, if fll.lSp holds for each compact interval [a, (3] in 

f±{x)=exp{2{g{x)±h{x))) (1.15) 

then ( 11. 9p holds with 

\Q{x) - g{x)\ < h{x) (1.16) 

Proof. By Theorem 11.11 for any positive continuous function, r]{x), on 
[a,P] supported on (a,/?), we have 

"^^^^ dx< [ r]{x)dfi{x) < [ -^P)-dx (1.17) 



Trf+{x) J J vr/_(x) 

from which absolute continuity of fi \ {a, [3) and fll.l4p are immediate. 
This in turn implies f ll.lSp and fll.l6p . □ 

Thus, we need to show a^p^ +Pn-i bounded as n — * oo, but with 
bounds that diverge 1 2. The difference equation is 

Pn+l \ I fx - bn+l -l\ ( Pn 



^A^+,{x)( ^" ] (1.18) 

Here 

det(A„) = 1 tr(A„) = x-bn (1.19) 

In like f 1 1.8 1) where 6„ is negative and monotone increasing, a 

fundamental object is the turning point, the integer, N{x), with 

x-bn>2 iin<N{x) (1.20) 
x-bn<2 if n>N{x) (1.21) 

If 7„(x) is defined by 7„ > and 

x-bn = 2 cosh(7„(x)) (n < N{x)) (1.22) 
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then one expects some kind of exponential growth as exp(^"^j^ •yj^x)), 
and we will prove that 

exp('^7,(a;)') < Pn{x) < {N + 1) exp(j2^,{x)] (1.23) 

As one expects, there is an intermediate region N{x) < n < Ni{x) 
and an oscillatory region n > Ni{x). We will see that so long as one 
is willing to accept 0((6Ar+2 — b^+i)^^) errors (and they will typically 
be very small compared to exp(^^j^ 7j(x))), one can actually take 
A^^i = + 2 (!) and use the method of proof for Theorem II .11 to control 
the region n > Ni. Thus, the key will be (11.231) and we will get (11.161) 
where 

N 

9{x) = J2^,{x) (1.24) 

i=i 

and 

hix) = 0(max(log(iV), log((6jv+2 - &7v+i)"'))) (1.25) 
The discussion of turning points sounds like WKB — and the reader 
might wonder if one can't obtain our result via standard WKB tech- 
niques. There is some literature on discrete WKB [6l [351 EZ] , but 
we have not seen how to apply them to this situation (for a different 
application to OPRL, see or, because of a double n ^ oo, x ^ 2 
limit, how to use the continuum WKB theory (on which there is much 
more extensive literature) to the continuum analog of our problem here. 
That said, the current paper should be regarded as a WKB-like anal- 
ysis. 

In Section [21 we discuss the case a„ = 1, 6„ < bn+i < 0. In Section [3l 
we discuss 6„ = 0, a„ < a„+i < 1. It is likely one could handle mixed 
dri"} cases with more effort. In Section [H we discuss some Schrodinger 
operators. Finally, in Section [3, we discuss examples including (11.81) 
and flTTT]) . 

It is a pleasure to thank Fritz Gesztesy, Uri Kaluzhny, and Doron 
Lubinsky for useful discussions. B. S. would like to thank Ehud de 
Shalit for the hospitality of the Einstein Institute of Mathematics at 
the Hebrew University where some of this work was done. Y. L. would 
like to thank Matthias Flach for the hospitality of the Department of 
Mathematics at Caltech where some of this work was done. 

2. Monotone 6„ 
In this section, we will prove: 
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Theorem 2.1. Let dji he the spectral measure associated with a Jacobi 
matrix having an = I and 

bn < bn+i <0 bn as n ^ oo 

Define N{x) for x in (0,2) and near 2 by fll.20p / (11.211) and 7„(a;) 
by (11.221) . Then dfi is purely absolutely continuous on (—2,2), where 
w = ^ is continuous and nonvanishing on (—2,2), 

Ci{x + 2) <w{x) <C2ix + 2)-^ for xe {-2,0] (2.1) 
and on (0, 2), 

w{x) = e-2'3(") (2.2) 

where 

\Q{x) - g{x)\ < h{x) (2.3) 

where 

N(x) 

^7(^) = E^^•(^) (2.4) 

and h{x) is given by 

e^(-) = CN{x){bNix)+2 - bNix)+ir\2 - x)5 (2.5) 
for an explicit constant C {dependent on sup|6„| but not on x). 

Remark. Typically, h is much smaller than g. For example, if 6„ is 
given by ([HD, g{x) = 0{{2-x)^~^) and e'^^^) = 0{N{xy+^{2-x)^) = 
0{{2 - a;)-(5+|))^ so h{x) = 0(log(2 - x)-^). 

As we explained in the introduction, we need to study the asymp- 
totics of Pn{x) as a; I 2 with some uniformity in n. Given that a„ = 1, 

Pn+i{x) = (e^"+^ + e-^"+')pn{x) - pn-i{x) (2.6) 

p^i{x)=0 Po{x) = l (2.7) 

which suggests we define for n < N{x), 

iJn{x) = e-^"=^^^Pn{x) (2.8) 

so ipn obeys 

V'„+i(x) = (1 + e-27n+i)^„ _ e-(7"+7n+i)^^_^ (2.9) 

V'-i(x)=0 tPf^{x) = l (2.10) 
Lemma 2.2. For < n < N{x), 

^n+l > i'n (2.11) 

In particular, 

M^) > 1 (2.12) 
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Proof. As a preliminary, we note that 6„ < bn+i imphes x — bn > 

X - 6„+i, so 

< 7n+l < In (2.13) 

By (ED, 

(2.14) 

For n = 0, - z^n-i = 1 > and Vn-i = > 0. By fl27[D and fl233D 
(which imphes e~'''"+i— e"'''" > 0), we see inductively that ipn+i—i'n > 0, 
and so, ipn+i > "^n > 0, proving (12.111) . □ 

Lemma 2.3. Define for n = 0,1,2, ... , N{x) — 1, 

Wn = e^-+'iljn-e-^-ijn-i (2.15) 

Then 

Wn < e^"+^ (2.16) 
Proof. Wq = e"'^ < e^^ , starting an inductive proof of fl2.16l) . By (12. 9p . 

so 

Wn+i = e(^"+2-^"+i)(Vr„ + e-^"+'^pn) - e"^"+iVn 

= e(^"+2~^"+i)W^„ + e-T"+i(e(^"+2-^"+i) - 1)^„ (2.17) 

< e^^"+^-^"+'^Wn (2.18) 

since fl213D implies e'^"+2 < e'^'^+i and ^„ > 0, (e'>"+2-T"-i - 1)^„ < 0. 
Thus, Wn < 6^"+'^ implies Wn+i < e'^"+'^ and (12.161) holds inductively. 

□ 

Lemma 2.4. For n = 0,1,2, ... , N{x) - 2, 

1pn+l<l+^n (2.19) 

So, in particular, for < n < N{x), 

^n<n + l (2.20) 

Proof By fl27[5D . 

since e-^"+W„+i < 1 by fl2J6|) and 7^ > implies e-('^"+i+'^"+2) < 1. 
This proves (I2.19p . which inductively implies (I2.20p . □ 

We summarize with: 
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Proposition 2.5. For any n with 1 <n < N{x), 

eE"=i7,W < <{n + l)eS?=i^^(^') (2.21) 

In particular, if 

Vn{x) = Pn-l{xf + Pn{xf (2.22) 

then 

e2E"=i7.W < ^^(^) < 2{n + ife^T.U'^ji-) (2.23) 

Proof. ((2211) is an immediate consequence of (ES]), fl2:T2l) and i^JM . 

□ 

Suppose X G (0, 2). For > N{x), define /€„(x) by < < | and 

a; - 6„ = 2cosft„(x) (2.24) 

so > hn+i > bn implies 

< K„ < Kn+l 

and bn — > implies 

ft^n ^ ft^oo = cos^^d) (2.25) 
For later reference, we note 

sin(«:oo) = (l-(f)')^ = |(4-x2)^ (2.26) 

So as X t 2, 

Koo = (2-x)^ +0((2-x)i) (2.27) 

We first present a matrix method following Kooman [12] to control 
the region [N{x) + 2, 00). At the end, we will discuss an alternate 
method using scalar Priifer-like variables. 

By ( ll.lSp . for n > N, An has eigenvalues e^*''". In fact, 

{'T V) (e-) (ei.) 

SO if 

yi^) = (X i) (2-29) 

and 

^ e''""^ 
then 

A„(x) = F(/€„)F(/s:„)r(/s:„)"' (2.31) 

Next, notice that 

Yi^r' = ( T-^^ ~]) (2.32) 

2t sm K V — e — i / 



V{K)=rn X) (2-30) 
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Following Kooman [12j, we write for n > £ > N{x), 

T„(x) = (2.33) 

= Y{K^)VnY{K„)-'Y{Kn^i)Vn-i ■ ■ ■ Y (Ke+l)-' 

and since ||Kt(/^)|| = 1, 

n-l 

\\Tn\\ < \\YM\\ m^e+.r'W n \\Y{^H^)-'Y{^,)\\ (2.34) 

j=e+i 

This prepares us for two critical estimates: 
Lemma 2.6. We have 

\\Y{^,^,r'Y{^,)\\ < 1 + . , \ ' (2.35) 
so, in particular, 

\\Y{k,^,)-'Y{k,)\\ < 1 + h±LZ^ (2.36) 
Proof. By (ESSl) and (I2J2]), 

^(S+i) ^(S) - 1 = 2sin(Kj+i) Ve-^''^ - e-^'^^+i e'''^ - e'''^+^ ) 

(2.37) 

If A = (aij) is a 2 X 2 matrix, 

\{ip,A^)\ < max(|ai,-|)(|(^i| + |<^2|)(|^i| + \H) 

< 2max(|a,,|)(|^in + \^2\')H\H' + \H't' 
since (|x| + \y\) < \/2{\x\'^ + \y\'^)K so 

\\YiK,+,)-'Y{K,) - 111 < . / |e-^+^ - e-^ l 

sm(Kj+i) 

which implies f l2.35p . 

(I2.35P implies (I2.36p since | > Kj+i > Hj implies sin(Kj+i) > sin^Kj). 

□ 

Remark. That (12.361) holds with a 1 in front of {kj+i — sin(Kj) 
is critical. Lest it seem a miracle of Kooman's method, we give an 
alternate calculation at the end of this section. 

Lemma 2.7. We have that 

n f 1 + ^^^±1^') < ^ exp{K^e{^^)) (2.38) 
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where 



e{y)= sup [-1—--] (2.39) 
o<x<y Vsin(a;) x/ 



Remark. Since sin(x) = x — \ + O(x^), = ~ + f + 0{x^) and 



6 ' ^"^ ' sin{x) X ' 6 

since sin(x) < x, we see e{y) is finite and 

e{y) = 0(1) as y i (2.40) 

Proof. We have 

^ <— + e(«oo) (2.41) 



so, since Kj+i > k 



'J' 



1 + < ^ + (/?,•+! - Kj)e(/s:oo) (2.42) 

sin(Kj) 

< ^(1 + («:,•+! -fi:,-)e(«:oo)) (2.43) 

< ^ exp((fi;j+i - /tj)e(Koo)) (2.44) 

from wliicli (12.381) is immediate if we note that Kqo — < fi^oo- D 

Proof of Theorem \2.1[ By (12.341) and Lemmas 12.61 and 12.71 if Tk{x) is 
the transfer matrix from N{x) + 2 to > N{x) + 2, then uniformly in 

||T„|| < 2(sin(/tAf(^)+2))"^ — exp(Kooe(/too)) (2.45) 

where we also used ||F(Kfc)|| < 2 and ||5^(/tAr(a;)+2)^"'"|| < 
2/2sin(K7v(^)+2)- 

As X T 2, ^ 0. Indeed, by ([22ZD, = (2 - x)^ + 0((2 - x)t). 
Moreover, by the definition of N{x), 

X - bN+i < 2 (2.46) 

while 

X - &Af+2 = 2 cos(k7v+2) (2.47) 

so 

2(1 - cos(KAr+2)) > bN+2 - bN+1 (2.48) 
Since A^(x) oo, bN{x)+2 ^ so /«Ar+2(x) — > and (12.481) implies 

^N+2{xf > (1 + o(l))(6^+2 - bN+i) (2.49) 

Thus in (|2.45|), [Kiv(^.)+2 sin(fi:jv+2)]~^ < (1 + o(l))(&7V+2 - Wi) and 
(I2.45P becomes 

sup ||f„|| < C{2 - x)-^{bN+2 - bN+i)-' = Aix) (2.50) 

n>N(x)+2 
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where now T„ transfers from — 1 to n and we use the boundedness 
from - 1 to + 2. Using 

(2.51) 

and (I2.23p . we obtain for all n > N, 

CiA(a;)-VS"^^(^') < {\pn\' + \Pn+if) < CA(a;)2Ar(x)2e2Ef7.W 

(2.52) 

which, given Corollary 11.31 implies fl2.2p - fl2.4p . 
In going from fl2.5ip to (12.520 . we used 

det(f„) = 1 ^ ||f-l = ||f„|| 

We also need to control the region x > —2 with 2 — x small. By 
replacing x by —x (and p„(x) by (— l)"p„(— x)), this is the same as 
looking at X + bn with still 6„ < bn+i < 0. We define 9n{x) by 

2cos(^„(a;)) =x + 6„ (2.53) 

so 

^1 > ^2 > ■ ■ ■ > ^oo = «:oo = (2 - s)5 + 0((2 - x)5) (2.54) 
As above, we have (I2.35p . so 

iim+irm)ii<i+%^ (2.55) 

but since Oj^i < 6j, we have 

l^i+i| ^j+i ^j+i 

and we find that, with T„ being the transfer matrix from 1 to n, 

This bound on the transfer matrix and Corollary 11.31 yield ( 12. ip . □ 

Remark. It might be surprising that (12. ip has (x + 2), (x + 2)^^ rather 
than (x + 2)5, (x + 2)~5 (because Carmona's bound relates w{x) to 
||T„|p and sup ||T„|| goes like (2 — x)2). Even in the free case, bounds 
from Carmona's formula give the wrong behavior: sin(n^) + sin^((n + 
1)^) have oscillations that cause the actual square root behavior in the 
free case, and bounds based only on ||T„|| lose that. 
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That completes the proof of Theorem 12.11 the main result of this 
paper. Here is an alternate approach to controlling for n > N, 
using the complex quantities: 

^n=Pn- e"*^"p„_i (2.58) 

so, since pj is real, 

sin(fi;„)|p„_i| = |Im(-$„)| 

< |$n| (2.59) 

By (EJlD, we have 

Pn+i = (e*''"+i + e-*''"+^)p„ - (2.60) 

so 

$n+i = e'''"+'[pn - e-*^"+^p„_i] 

= e'^^+'^n + e'^'^+'ie-'''" - e-'^"+')pn-i (2.61) 

Using (12391) . 



I<fn+l|<|$n| + ^^.7^|$n| (2.62) 

sm U« 



and similarly. 



I<^'n+l|>|$n|-^^^^r7^|<fn| (2.63) 

These replace (I2.36p and imply, via Lemma 12.71 and the analysis in 
([236]), that 

Ci{2 - x)-kbN+2 - bN+i) < < C{2 - x)-^{bN+2 - bN+ir' 

\^N+2\ 

Since 
and 

2|$„|2 > sin2(/€„+i)(|p„p + bn-iH 
we can go from this to Theorem 12. 1[ 

3. Monotone a„ 
In this section, we will consider 

bn = ctn+i < a„ < 1 a„ ^ 1 (3.1) 

The weight will be symmetric, the measure purely absolutely continu- 
ous (i.e., no eigenvalues outside [—2,2]), and so for non-Szego weights, 
the integral will diverge at both ends. Here is the main result: 
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Theorem 3.1. Let dfi{x) = w{x) dx be the measure associated with 
Jacobi parameters obeying fl3.ip . For any x G (—2, 2), define N{x) by 

2an < |x| forn<N{x) 2an > \x\ forn>N{x) (3.2) 

and 7n(a;) for n < N{x) by 

^ = 2cosh(7„(x)) (3.3) 
an 

Then 

w{x) = e-2Q(") (3.4) 

where 

\Q{x) - g{x)\ < h{x) 

N{x) 

and h{x) is given by 

gMx) = CN{x){aN(.)+2 - aNi.)+i)-' (3.5) 



The proof wiU closely mimic the proof of Theorem 12. ![ so we will 
only indicate the changes. By symmetry, without loss, we can suppose 
a; > 0. The recursion relation becomes 

Pn+i{x) = (e^"+i(^) + e-^"+^(^))p„(x) - ^p„_i(x) (3.6) 

where we note, by (13.31) . that 

Qn ^ cosh(7„+i(a;)) 
a„+i cosh(7„(a;)) 
Define ipnix) by (12. Sp . so (12. 9p becomes 

V„+i(x) = (1 + e-27"+^(^'))^„(x) - ^ e-(^"(^)+^"+^(^)Vn-i(x) (3.8) 
f l2:T0D still holds. 

Lemma 3.2. ■i/'ri+i > V'n; so ipni^) > 1 /or n > 0. 
Proo/. We still have fl2J3D . and flCTj) becomes 

(3.9) 

Since a„ < a„+i, < 1, and so 
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Thus, by (12. 9p . ipn+i — V^n > and ipn+i > inductively. □ 
Lemma 3.3. 

e-- < e-- (3.10) 
cosh(7„+i) 

Proof. This is equivalent to 

g7n + 2+7n + l _j_ g7n + 2-7n + l <; g7n + 2+7n + l _j_ g77i + l-7n + 2 ]^]^^ 

SO to 7„+2 - 7n+i < 0, SO to (I2.13p . □ 
Lemma 3.4. Define 

Wn = e^"+'^Pn - e-^-^n~x (3.12) 

< e^"+i (3.13) 

Proof. (I3.13P holds for n = by (I3.12p for = 0, so we can try an 
inductive proof. The analog of (12.170 is 

Wn+i = e(T"+2-^"+i)vr„ + e-^"+i ( e(^"+2-'^"+^) - ^ ) ^„ (3.14) 

\ O'n+2 J 

By (I32D and ([SlOD, 

g{7n+2-7n + l) _ ^"+1 < Q 

SO (I3.14P says 

by induction. □ 
Lemma 3.5. ipn+i < 1 + "^n so inductively, ipn ^ n + 1. 
Proof By ( KWl and fim]) . 

O'n+2 

since ^ < 1. □ 

O.n + 2 — 

If now 

Vn{x) = Pn-l{xY + alpn{xf (3.15) 

then we have proven (I2.23P for large n. 

To control the region n > N{x) + 2, we use the scalar variable 
technique from the end of Section [21 Define k„ for n > N{x) + 1 by 
(recall x > 0) 

X 

— = 2cos(k„(x)) (3.16) 
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SO an+i > Qn impHes 



Define 
Then 

Lemma 3.6. (i) 



\Pn-l\ < - 



1$. 



n 



I'^n+ll 
l$nl 



< 1 + 



<1 + T 



le**^" cos{K,n) - e"^"+i cos(k„+i) 



COS[K,n) Sm[K,n 



Proof, (i) This comes from |Im$„| = sin(/t„)(p„_i) 
fii) From 



we obtain 



By (EIED, 



|$„+i -e^'""+^$„| 



■Pn-l 



Pn-l 



a„ COS(K„+iJ 



(3.17) 
(3.18) 

(3.19) 

(3.20) 
(3.21) 



(3.22) 



(3.23) 



a„+i cos(k„) 

so (IX^ and (IXT^ imply ([220]) • This in turn implies flS:^ since 

e'""- cos(/t„) - e^'^'^+i cos(fi:„+i) = \ {e^'"" - e^'^"+') (3.24) 

□ 

With this formula, we can mimic the proof of Theorem 12.11 to com- 
plete the proof of Theorem 13. 1[ 

4. SCHRODINGER OPERATORS 

In this section, we consider Schrodinger operators H = — ^ + V{x) 
on L^([0, cxd)) where one places u{0) = boundary conditions. H is 
unitarily equivalent to multiplication by E on L^(M, dfi{E)), where 
is the conventional spectral measure (see [Sj [IHl [23] ) . If u{x,E) obeys 

~u" + Vu = Eu u{0,E) = 0, u'{0,E) = 1 (4.1) 
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then Carmona's formula [1] takes the form 

IT fii J uj , , , 

In particular, if uniformly in compact subsets of E E (0, oo), 
exp(2(^(E) - h{E))) < liminf + 

x— »oo 

< limsup {\u{x,E)\'^ + \u{x,E)\'^) 

< exp{2{g{E) + h{E))) (4.3) 

then dfx is purely absolutely continuous on (0, oo), dn(E) = e"^*^*^^) dE, 
and 

\Q{E) - g{E)\ < h{E) (4.4) 

We want to assume the following conditions on V: 

(a) V is Ci on [0,oo). 

(b) V is positive and strictly monotone decreasing on [0, oo). Indeed, 

V'{x) < (4.5) 

(c) 

lim V{x) = (4.6) 

X— >oo 

Of course, the canonical example is 

V{x) = {x + xoy^ (4.7) 

Our main result in this section is: 

Theorem 4.1. Let V obey (a), (b), (c) so dn{E) = e'^^^^^ dE. Define 
forE < V{0), 

N{E) = V-\E) 



(4.8) 



so 

V{x) > E if x< N{E) 
V{x) <E if x> N{E) 
For X < N{E), define 

j{x,E) = {V{x) - E)^ (4.9) 

Then holds where for E <V{0), 

rN{E) 

g{E) = / -f{x,E)dx (4.10) 



and for E <V{0), 

^h(E) ^ CN{E) {V{N{E)) - V{N{E) + l))-^^^ (4.II) 
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This proof will illuminate the proofs of the previous two sections. 
We begin with an analysis of the region x < N{E). We define 

4.12) 

4.13) 

4.14) 
4.15) 

then 

□ 

4.16) 
4.17) 
4.18) 



ipix) =u{x,E)expy- J -f{y,E)dy 
and are heading towards 

< ^'{x) < 1 

Lemma 4.2. For < E < V{0) and x < N{E), we have 

(a) u{x) > 1 

(b) u{x) > X 

Proof, u" = '-f'^u, so u" > 0. This implies u'{x) > u'{0) = 1, and 
Lemma 4.3. For E < V{Qi) and x < N{E), 

i)\x) > 

Proof. Let 

/(x) = u'{x) — 7(x)u(x) 

so 

iP'{x) = f{x)exp(^- j{y,E)dy 

and (I4.16P is equivalent to / > 0. Note that 

f + if = ^" ~ 7^' ~ I'u + lu' — 
= —j'u 

since (14. ip says 

u" = 

(14. 5p implies 

i'{y) < 

so (^J^ says 

which, given /(O) = 1, implies / > and so tp' > 0. 
Lemma 4.4. Let 

W{x) = ip'ix) + 2-f{x)ilj{x) 
Then W'{x) < and so 

ij'{x) < 1 



4.19) 

4.20) 
4.21) 

4.22) 

□ 

4.23) 
4.24) 
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Proof. By (Km . 



so 



W\x) = {u" + 7m' + 7'm - 7m' - -f\)e- ^^v) <^y 
= 7'Me" -fo 'y(y) '^y 
< 

by Km . But W{x = 0) = ^'(0) = 1, so 

W{x) < 1 

and thus 



(4.25) 



(4.26) 



(4.27) 

(4.28) 

□ 



Proposition 4.5. If E is such that N{E) > 1, then 

g-2y(o)g2 < uiNiE)y+u'iNiE)y < (iV(E)2+l)e2/o"^"'^(2^)'^J' 

(4.29) 

Proof. Since 7(A^(^)) = 0, 

tp'iN{E)) = u\N{E))e'-^^^^''^^^y'^'^y 

so Q <ip' <1 and V^(0) = yield tlie upper bound in (14.291) . 

For the lower bound, fl4.15p implies u{l) > 1. So, since '~f{y) < 
7(0) < V{0), 

V'(l) > e-^(°) (4.30) 
which, given that ip' > and N{E) > 1, implies 

«(iV(E))>e-^(°)e/o"'"^(J')'^^ □ 

In the region [N{E), N{E) + 1], we note that since 



1 V{x) - E 
1 



< 1 + \E\ + \V{0)\ 



the matrix form of the Schrodinger equation implies that if C{x) = 
|m(x)P + |u'(x)p, then 

e-^(^+\^\+nm-y\c{y) < C{x) < e2{i+l^^l+v'(o))|:.-s/lc'(^) 

giving a constant term in e'^^^^ in (14. lip . 

Finally, in the region [N{E) + 1, 00), we use the method of Appendix 2 
of Simon [30j (see also Hinton-Shaw [9]). Define for x > N{E), 

k{x,E) = ^yE-V{x) (4.31) 
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and define 



If 



u±{x) = expl ±.i / K{y) dy 
V Jn{e) 



F{x) ='-V'{x){E-V{x))-^ 



(4.32) 



(4.33) 



and if a{x), b{x) are defined by 

u{x) = a{x)u+{x) + b{x)u^{x) (4.34) 
u'{x) = a{x)u'_^_{x) + b{x)u'_{x) (4.35) 

tlien u" = —H?u is equivalent to (see Problem 98 on p. 395 of [28j ) 



a{x)\ ' 
b{x 



M{x) 



a{x] 
b{x) 



wliere 



witli 



M{x) = w{x 



,-1 



-F{x) ul{x)F{x) 
uI{x)fIx) -F{x) 



w{x) = u_^^{x)u^{x) — u'_{x)u^{x) 
= 2m(x) 

Proposition 4.6. Let M{x) be given by ( I4.37p . Then 

k(oo, E) 



\\M{x)\\ dx < log 



'N{E)+1 

Proof. Since = 1, 



k{N{E) + 1,E) 



\\M{x)\\ < \w{x)\-^ 



\F{x)\ \F{x)\ 
\F{x)\ \F{x)\ 

2\w{x)\-^\F{x)\ 



'V'ix)iE-Vix))-^ 
j-\og{{E-V{x))^ 



from which (14.391) follows. 
Proof of Theorem \4.1\ Let 

Y{x) 



u^{x) u^{x) 
u'j^ {x) u'_ (x) 



(4.36) 



(4.37) 



(4.38) 



(4.39) 



(4.40) 

□ 

(4.41) 
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Let T{x,y) be the (^,) transfer matrix from x to y and T{x,y) be the 
(2) transfer matrix. For y > N{E) + 1, we have just seen 



||T(iV(i?) + l,y)|| <exp / \\M{x)\\dx 

\Jn{e)+i 

K,{oo, E) 



oo 



(4.42) 



KiN{E) + l,E) 
On the other hand, 

\\Y{y)\\ <1 + K<2 (4.43) 

for K small while 

\\Y{yr'\\ = \det{Y)-'\ \\Y\\ < ^{y)-' (4.44) 

and 

T{x,y) = Y{y)f{x,y)Y{x)-' 

so 

Since E = V{N{E)), 

k{N{E) + 1, Ef = V{N{E)) - V{N{E) + 1) (4.46) 

and we have the bound (14.41) with the error built from e~^'^'^\ N(E), 
(K39\f . and (K^ . □ 

It is interesting that the differential equation methods of this section 
lead to terms that are identical to what we found in the discrete case. 

5. Examples 
We start with the continuum case. 
Example 5.1. 

V{x) = Cox-f^ l3<2 (5.1) 
Technically this does not fit into Theorem 14.11 since ^(0) = oo, but 
when /3 < 2, it is easy to extend the analysis. The spectral measure is 
g-2Q(£) ^j^g^g holds. 

mE)=(^) ' (5.2) 



V{N{E)) - V{N{E) + 1) ~ V'{N{E)) 

~ N{E)-^V{N{E)) 

= EN{E)-^ (5.3) 
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SO h{E) = 0{\og{N{EfE'^)) = 0{\og{E)). On the other hand, let- 
ting y = x/N{E), 

l-N{E) 

g(E) = / {V{x)-E)idx (5.4) 



= N{E)E^ / {y-^-l)^dy (5.5) 

1 , 

= E^N{E)p-^ {l-u)^u0~^du 
Jo 

Y(^)T(- - -) 
= e'^N{E)P'' (5-6) 

^ ^ r(i + i) ^ ^ 

using slU = y^ change of variables. Thus, 

1 rfSirfi _ l\ 

g(E)=c^^Ei-i c. = r' r(l + l) <5-') 

Since /? < 2, g{E) — > oo and is much larger than the log(i?) error. 
/3 = 1, the Coulomb case, has g{E) = CoCiE~2 and P = \, the quasi- 
Szego borderline, has g{E) = C^ciE"^. We emphasize that g occurs 
in an exponential, so w is very small near E = 0. □ 

Example 5.2. 

V{x) = Co{x + Xo)-'' (3<2 (5.8) 

We claim that the changes from Example 15.11 are small compared to 
log(i?) errors in h\ explicitly, 

g{E) = ciC^E^^-^ +0{1) + 0{e'^) (5.9) 

For in this case, 

N{E)= ' -xo (5.10) 



Co, 

and one changes variables toy = {x+xo)/{N{E)+xo), so (15. 5p becomes 
g{E) = N{E)e"^ [ {y~P - l)h dy (5.11) 

Js{E) 

where 
Then 

s{E) 



N{E)e'2 {y-^ -1)^ = N{E)E'W{s{Ey- 

Jo 
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= 0(1) (5.13) 

by fISJOD and (KT2^ . so 

g{E) = ciN{E)Ei + 0(1) 

= CiCo^E5-i +0(1) + 0(£;5) (5.14) 
as claimed. □ 

Now we turn to the discrete case. 
Example 5.3 (= (iLSi) ). 

a„ = 1 bn = -Cn~<^ (5.15) 



Define 



X 



5„ = Cn-^ - 6 (5.16) 



so 



x-b^ = 2 + 6n (5.17) 
We have (with [y] = maximal integer < y) 

N{x) = [{C-^6y^] (5.18) 

We have 6;v+2 - &jv+i = O(iV-^-i), so the RHS of is of order 
CiV(a;)^+^5i = 0(5^2-^) and thus, h{x) = 0(log(2 -x)) and we need 

to compute g{x) = Xl^^i'' 7i(^) 0{\og6) terms. 

We will suppose below that C < 1 and explain at the end what to 
change if C > 1. 

Define q to be the Taylor coefficients in 

oo 

cosh-^(l + I) = (5.19) 

e=o 

so, courtesy of Mathematica, 

1 3 5 

Co = 1 Ci = -— C2 = — -- C3 



24 640 7168 

and, for example, 

C20 = 34, 461, 632, 205/12, 391, 489, 651, 049, 749, 040, 738, 304 

(assuming that we managed to copy it without a typo). Thus, 

oo N{x) 
1=0 j=l 

Notice that since S > 0, 

Sj < Cj-^ (5.21) 
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SO, if + |) > 1, a crude (5-independent bound of Yl!j=i ^^^'^ can 
be summed independently of N{x). Moreover, if F is the function in 
( Km . then 

2^— = , 5.22 

so the Q power series has radius of convergence 4 and so < oo. 

Thus, if 

^0 = - i] + 1 (5.23) 

then 

oo N{x) ^ oo \ oo 

ENE^r^^lENjEr^^'""'^ (5.24) 

l=lQ j = l ^0 ^ j = l 

(since C < 1) so 

N N 

E^^= E Q 5^ 5^^+0(1) (5.25) 

If £ = — I occurs, then 

^ 1 1^1 ^ ^ 



AT 



TV 



O(logA^) (5.26) 



On the other hand, if £ < ^ — ^, then 

N N 

i=i i=i 

N 

= c^+i5^r^(^+^Hi-(imUo(i) 
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N 



•1 



(5.27) 
(5.28) 

C^+^/3"i7Vi-(^+i)/3 / ' «(i-^-f)(i - uY^-^ du + 0(1) 

(5.29) 



' -^iVM^+|)/^ + 0(l) 



In the above, fl5.27p comes from the fact that the function in the 
integrand is monotone decreasing, and if f{x) is monotone, then 



m > 



f{y) dy>f{j + l) 



so 



and 



N-l 



f{y)dy>Y,fU) 

j=2 



N 



</(l) 



(5.30) 



f{y)dy-J2f{j) 
i=i 

f l5.28p is the change of variables u = (ff)^- Finally, (15.291) comes from 
the same cancellation that occurred in fl5.13p . 

Since \N - ch'^ < 1 and < 1 - (£ + i)/? < 1, 



(5.31) 



Thus, we find 



c—^-ii—i — ir^+^+^ + o(iog5) 



o<^<(^-f) 



r(i + i) 



(5.32) 

If C > 1, we should not expand the power series of cosh"^ for small 
j (actually, as noted, the power series has radius of convergence 4 so 
we need only worry if C > 4). Instead, we do not expand for those 
j with Cj~^ > 1. That is only finitely many terms, so it adds 0(1) 
errors to We add back these small j terms to (15.250 . again 

making 0(1) errors. The final result does not change. □ 
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Finally, we will explore examples that lead to Q's roughly of the type 
f ll.lip to link to work of Levin-Lubinsky [IB] • We suppose 

a„ = 1 - /(log(n + 1)) (5.33) 

where the /'s we have in mind are typically 

/(x) = (l + x)-" (5.34) 

or 

/(x) = logfc(a; + Cfc) (5.35) 

an iterated log (where Ck is chosen to keep all log's that enter positive). 
We will need 

Proposition 5.4. Let f be defined and C"^ on [log2,oo) and obey 



(i) 


fix) > 0, fix) < 0, fix) > 


(5.36) 


(ii) 


lim fin) = 

n— »oo 


(5.37) 


(iii) 


lim N%-f'i\ogN))^ = oo 

A''--»oo 


(5.38) 


(iv) 


r A- -f'{{^-e)k)\ ^ 

lim lim sup — = 1 

^10 V fe-oo -f'ik) J 


(5.39) 



Let 



N 

Sn = Y1 v//(logj)-/(logiV) (5.40) 
i=2 

Then 

lim ^ , = ^ (5.41) 

N^oo Ni-f'i\ogN))2 2 

Remark. It is easy to see that if fix) = e~^^ (i.e., /(log(?7, + 1)) ~ 
in + 1)-'=), then flCTll fails. In this case, both (15:381) and fICTD fail, 
but they hold for the /'s of (ESI and fl5^ . 

Proof. Since (—/')' < and if a; < 



/(a;)-/(y)= / (5.42) 

we have, 

(2/ - x)i-f'iy)) < fix) - /(y) <iy- x)(-/'(x)) (5.43) 
We thus get a lower bound 

/(logj) - /(logiV) > (-/(logiV))(-log(i)) (5.44) 
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SO 

N 

Sn > iV(-/'(logiV))i Yl F(-log(if))^ (5.45) 

i=2 

As oo, the sum converges to /q^(— log(a;))5 dx = ^ (courtesy of 
Mathematica) . Thus, 



hmmf(LHS of dEUD) > ^ (5.46) 

For the upper bound, fix £ > and break Sn = + S^^^ where 
S^j^^ has 3 < N^-' and ^JJ^ has j > N^-' . Clearly, 

4')</(log2)iVi"^ (5.47) 

so, by hypothesis f l5.38p . it contributes to the ratio in (15.411) as — >• 
oo. 

For we use the upper bound when j > N ^ 

/(logj) - /(log AT) < -/'((I - e) \ogN){- log(i)) 
which yields (since the Riemann sum still converges to the integral) 



limsup(LHS of ^M) < — limsup ' ^'^^^ ^^^^ 



2 fc— >oo V /'(^) 

Since e is arbitrary, we can use (15.391) to complete the proof of fl5.4ip . 

□ 

Example 5.5. Let a„ have the form fl5.3ip where / obeys all the 
hypotheses of Proposition [531 By (13. 2p and (13. 3p . A^(x) roughly solves 

^ " (5.48) 



l_/(log(iV + l)) 
namely, 

Ar(x) = [exp(ri(l- !))]-! (5.49) 
For example, if / is (I5.34p . then 

N{x) = [exp((l - f)-" -!)]-! (5.50) 
Next, define z by ^ = 1 + |, namely, 

z = ^-2 (5.51) 
where f > 2. Writing x = 2 — 6 and a = 1 — /, we see 

z = -6 + 2f + Oif) + Oif5) (5.52) 
Taking into account that N{x) is such that 

2/(log(Ar + 2))<5<2/(log(Ar + l)) 
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and that (15 .191) says 

cosh-i(l^) = v^ + 0(2t) 

we see that 

7,(x) = v/2/(log(j + l))-5 + 0(/i) + 0(/^5) 

and thus 

N(x) 

i=i 

is asymptotically the same as \/2Sn- Thus, 

\Q{x) - g{x)\ < h{x) (5.53) 

where 

g{x) = y|iV(x)(-/'(logiV(a;)))^(l + o(l)) (5.54) 

and 

h{x) = 0{\ogN{x)) + 0{\og{l-l)) 

N{x) is huge, so while log A^(a;) ~ (1 — f)~" in case (15.341) . it is still 
small relative to g{x). 

The reader may be puzzled in comparing our results with those of 
Levin-Lubinsky [18] . They have no and their relations (after mak- 
ing the modifications from [—1, 1] to [—2,2]) suggest 

l-a„ = (logn)-^(l + o(l)) (5.55) 

should correspond to 

Q{x) = exp((l - f )-") (5.56) 

so there is no sign of (-/'(logA^(x)))2 either. 

The mystery is solved by the fact that multiple Q's lead to the same 
leading asymptotics for a„. In their scheme, after corrections to move 
to [—2, 2], leading asymptotics for / are given by 

n = Q(l-2(/(n)(l + o(l)))) (5.57) 

If 

Q(x) = ei/(^-t) (5.58) 

then 

n = exp((/(n))-^) (5.59) 

solved by 

/(n) = j^ (1 + 0(1)) (5.60) 

Changing fl5.58p to 

g(x) = f(l-f)exp((l-|)-i) 
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is solved by 

f{n) = l/(log(f logn) + 0(loglogn)) 

Since 

log ^ logra = logra + log2 n + log(f ) 
fl5:60|) still holds! □ 
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